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Abstract 

We analyze the effects of finite-range corrections in halo effective field theory 
for S-wave proton halo nuclei. We calculate the charge radius to next-to- 
leading order and the astrophysical S-factor for low-energy proton capture 
to fifth order in the low-energy expansion. As an application, we confront 
our results with experimental data for the S-factor for proton capture on 
Oxygen-16 into the excited 1/2^ state of Fluorine-17. Our low-energy theory 
is characterized by a systematic low-energy expansion, which can be used 
to quantify an energy-dependent model error to be utilized in data fitting. 
Finally, we show that the existence of proton halos is suppressed by the need 
for two fine tunings in the underlying theory. 

Keywords: halo nuclei, charge radius, radiative capture, effective field 
theory 


1. Introduction 

The quantitative description of both nuclear structure and reactions on 
the same footing is a major challenge of contemporary nuclear theory. With 
new experimental facilities such as FRIB and FAIR at the horizon, the task to 
find improved approaches for nuclear reactions has become very urgent. Ab 
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initio approaches to calculate nuclear scattering observables are limited by 
the computational complexity of the nuclear many-body problem. Scattering 
models perform well but use a number of uncontrolled approximations that 
make the errors of such calculations difficult to quantify. 

Faced with these problems, it is important to note that there are a number 
of systems in the chart of nuclei for which the effective number of degrees- 
of-freedom is signihcantly smaller than the number of nucleons. This phe¬ 
nomenon is known as clustering, with alpha clustering, e.g. in the Hoyle 
state of being the most prominent example. Clustering becomes even 
more extreme for so-called halo nuclei which consist of a tightly-bound core 
nucleus and a few weakly-bound valence nucleons [1, 2, 3, 4]. This reduc¬ 
tion in the number of degrees of freedom is the signature of a separation of 
scales in the system. In the case of a one-nucleon halo nucleus, the scale sep¬ 
aration is manifest in the small ratio of the one-nucleon separation energy 
and the binding and excitation energies of the core. For typical momenta 
on the order of the one-nucleon separation energy, it allows for a systematic 
low-energy expansion in the ratio of these two scales. This expansion can 
then be employed to calculate nuclear observables in a model-independent 
and systematically improvable manner. This approach is called halo effec¬ 
tive held theory (Halo EFT) [5] when a held-theoretical approach is used for 
the construction of the interaction and the calculation of observables. Halo 
EFT employs the minimal number of degrees-of-freedom (core and valence 
nucleons) and parameterizes the interaction in terms of a few measurable pa¬ 
rameters. In addition, so-called core polarization effects become important 
if the core has low-lying excited states. These can be taken into account 
by including excited states of the core as explicit degrees of freedom in the 
effective theory. 

Neutron halo nuclei occur rather frequently in the chart of nuclei along 
the neutron dripline and have been studied in Halo EFT [5, 6, 7, 8, 9, 10]. 
Proton halo systems exist too, but are less common due to the delicate inter¬ 
play between attraction from the strong interaction and the repulsion from 
the Coulomb interaction. The effects of the Coulomb interaction were hrst 
included into an EFT with contact interactions by Kong and Ravndal [11]. 
Proton halos were considered recently in Refs. [12, 13, 14]. The Coulomb 
interaction introduces a new scale into the problem that can be understood 
as a result of the presence of a Coulomb barrier. This new scale is refered to 
as the Coulomb momentum and it is given by the inverse of the Bohr radius 
of the system. The introduction of a Coulomb momentum can complicate 
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the power counting since it interferes with the separation of scales. Higa et 
al. [15], e.g., treated the Coulomb momentum as a high-momentum scale in 
their study of a —a scattering. However, this treatment is not always appro¬ 
priate. The correct scaling of the Coulomb momentum will always depend 
on the system to be considered. 

In this paper, we will extend the calculation performed in Ref. [12] by in¬ 
cluding higher-order effects due to the hnite range of the interaction between 
core and proton. We also consider higher-order electromagnetic interactions. 
Specihcally, we will consider the charge radius of S-wave halo nuclei and ra¬ 
diative proton captnre into a halo state. Our analysis of hnite-range effects 
also addresses the qnestion of why there are more nentron halos than proton 
halos in nature. 

This manuscript is organized as follows. In Sec. 2, we introdnce Halo EFT 
for S-wave systems and explain how the Coulomb interaction is inclnded into 
calculations. In Sections 3 and 4, we apply Halo EFT to calculate the charge 
radius of proton halo nuclei and radiative proton captnre, respectively. We 
address the aforementioned issne of hne tnning in proton halo nnclei in Sec. 5. 
Results for the excited 1/2’'' state of Fluorine-17 are presented in Sec. 6. In 
particular, we extract the threshold S-factor for radiative proton capture on 
into ^^F* and the corresponding asymptotic normalization coefficient. 
To this aim we employ an order-by-order fit to experimental radiative cap¬ 
tnre data at finite energies and we demonstrate how to qnantify theoretical 
nncertainties within Halo EFT. Finally, we summarize our findings in Sec. 7. 


2. Theory 


In Halo nnclei, the core and the valence nucleons are the effective de¬ 
grees of freedom. The Halo EFT Lagrangian can therefore be constructed 
using only a core and a nncleon field. For proton halos this was first done 
in Ref. [12]. In order to simplify the inclusion of finite range effects and 
fntnre extensions to two-proton halos, we use an eqnivalent approach that 
introdnces an auxiliary halo field d with the quantum numbers of the halo 
nucleus [5], leading to the Lagrangian 
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Figure 1: The full halo propagator defined iteratively. The full halo propagator is denoted 
by the thick double line. The thin double line denotes the bare halo propagator, the 
dashed single line denotes the core held, the solid single line denotes the proton held and 
the shaded blob the Coulomb four-point function y dehned in Fig. 2. 

The Lagrangian including nucleon and core fields only can be obtained by 
integrating out the halo field using the classical equations of motion. In 
Eq. (1), ipQ denotes the proton field with mass mo and ipi the core field with 
mass mi, while g and A are the low-energy constants of the theory and 
^tot = + ^ 2 - The parameter z/ = ±1 allows for the effective range to 

be both negative and positive, since we define the coupling g as purely real.^ 
The ellipsis denote operators with more fields and/or additional derivatives. 
The covariant derivative is = <9^ -l-ieQA^, where Q is the charge operator 
and e > 0 is the elementary electric charge. Low-lying excitations of the 
core can be included explicitly in the calculation by introducing additional 
core fields (See, e.g.. Refs. [7, 13, 14] for more details). In the calculations 
below, we do not need to keep track of the proton spin since the core-proton 
interaction does not change the spin of the proton. 

It is important to assign proper scaling dimensions to the fields and oper¬ 
ators of the Lagrangian (1) such that we can organize the theory in a power 
counting. The proton and core fields have the scaling dimension ~ 3/2 
and the scaling dimension of the halo field is d ~ 2 (see for example [7] for 
a discussion of the scaling dimension of such auxiliary fields). The covariant 
derivatives scale as Dq ~ 2 and Dj ~ 1. Note that we do not explicitly 
count powers of mass factors since mass is not equivalent to energy in non- 
relativistic physics. 

The Feynman rules for the interactions of the Lagrangian (1) are as fol¬ 
lows: 


(i) The vertex factor for the strong contact interaction is —ig. 


^Note that the kinetic term of the halo field d has the wrong sign for ly = —1 and 
becomes a ghost field. This introduces no pathologies here since the d field never appears 
in loops. 
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(ii) The vertex factors for the Aq interaction with the proton, core, and 
halo held are —ie, —ieZc and —iue^Zc + 1), respectively, where Zc is 
the proton number of the core. 

(iii) The vertex factors for the interaction of the vector photon, Ai, with 
a proton, core, and halo held carrying incoming momentum p are 
—iep/nio, —ieZcp/mi and —iue^Zc + l)p/Mtot, respectively. 


From the Lagrangian (1), the one-particle propagator can be deduced to 


be 


iSk{po,p) = i 


Po - 


2mk 


%E 


-1 


( 2 ) 


For convenience, we will also dehne the proton-core two-particle propagator 


*5'tot(po,p) 



P^ 


2m 


+ ie 


R 


-1 


( 3 ) 


where ttir denotes the reduced mass of the proton-core system. In this work, 
we will only need the halo propagator for a halo held at rest: 




i 

A + u {E + ie) 


( 4 ) 


The corresponding propagator for hnite momentum p can always be obtained 
by replacing E ^ E — p^ j (2Mtot)- The power counting for systems interact¬ 
ing through a large scattering length requires that the S-wave interaction is 
summed up to all orders [16, 17, 18]. The resulting full halo propagator is 
thus given by the integral equation shown in Fig. 1. For a halo held at rest, 
we obtain 


iD{E) 


i 

A + V {E + ie) + S(i7) 


( 5 ) 


The irreducible self-energy, S, includes strong and Coulomb interactions and 
will be discussed below. 

We include the Coulomb interaction between the core and the valence 
proton through the full Coulomb Green’s function. 


(k|Gc(i7)|p) = -^tot(i?,k)x(k,p;i7)5tot(i?,p) , (6) 

where p and k are the relative incoming and outgoing momenta and E is the 
energy. Here, y is the momentum-space Coulomb four-point function in the 
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Figure 2: The four-point function x defined iteratively. The wiggly line denotes a Coulomb 
photon exchange. External propagators are amputated. Otherwise, the notation is as in 
Fig. 1. 


center-of-mass frame of the proton and the core, dehned recursively in Fig. 2. 
To distinguish coordinate-space from momentum-space states we will denote 
the former with round brackets, i.e. |r). The Coulomb Green’s function can 
be expressed via its spectral representation in coordinate space 


(r|Gc(i!)|r') = j 


d^p 

(27r)^ E — p2/(2771^) -b ie ’ 


(7) 


where we dehne the Coulomb wave function through its partial wave expan¬ 
sion 

= £(2^ +■ r) . (8) 

1=0 7 * 

Here we have introduced p = pr and the Sommerfeld parameter p = kc/p, 
with the Coulomb momentum kc = ZcarriB., where a is the hne structure 
constant. We have also introduced the pure Coulomb phase shift ci = 
aigT{I + 1 + ip). For the Coulomb functions Fi and Gi, we use the con¬ 
ventions of Ref. [19]. The regular Coulomb function Fi can be expressed in 
terms of the Whittaker M-function according to 




(9) 


with the Ai dehned as 


Mv) = 


|F(/ + 1 + ip)\ exp [—7r?7/2 — i{l -b l)vr/2] 


( 10 ) 


2(21 + 1 )! 

We will also need the irregular Coulomb wave function, Gi, which is given 
by 

Gi{p, p) = iFi{p, p) + Bi{p)Win,i+i/2{2ip) , (11) 

where W is the Whittaker W-function and the coefficient Bi is dehned as 

exp {np/2 + iln/2) 


B,{v) = 


arg F(/ + 1 + ip) 


( 12 ) 



Figure 3: Feynman diagram for the irreducible self energy. The injected four-momentum 
is (i?, 0). External legs are amputated. Otherwise, the notation is as in Fig. 1. 


It is important to note that since we consider both free and bound states the 
absolute value and the argument of the F-function are given by 

|r(/ + l + ir])\ = A/r(/ + 1 + ir])T{l + 1 - ir]) (13) 


and 

+ i + + !!!,)! ■ (14) 

To obtain the fully-dressed two-particle propagator D, which includes 
strong and Coulomb interactions, we calculate the irreducible self-energy 
shown in Fig. 3. Using Eq. (6), it can be expressed by the momentum-space 
integral 

/ U3iL U3iL 

{k,\Gc{E)\k,) , (15) 

which can be written in coordinate space using Fourier transformations: 


lEiE) = - ig\0\Gc{Em 

= _ • 2 /“ d^p ^p(o)^;(o) 

J (27r)3 E - p2/(2mR)-b ie ‘ 


We evaluate this integral using dimensional regularization in the power di¬ 
vergence subtraction (PDS) scheme [17] and the result is [11] 

S(E) = , (17) 

TT 

with 

1 

H{ri) = + — - log {ig) , (18) 

Zig 

where ipiz) = T’{z)/T{z) is the logarithmic derivative of the Gamma func¬ 
tion. The divergent part is given by 


j.div ^ 


g'^kcm^ 
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+ log 


2kc I 


+ 1 - 




+ 


9 mn/j, 

271 


(19) 


7 



Figure 4: Scattering amplitude for elastic proton-core scattering. The notation is as in 
Fig. 1. 


with the renormalization scale /i. Note that is independent of energy and 
therefore will vanish when we take the energy derivative of the irreducible 
self-energy to arrive at the LSZ residue below. 


2.1. Renormalization 

Expressions for EFT low-energy constants such as g and A dehned above 
are frequently obtained by matching them to elastic scattering observables. 
The amplitude for elastic proton-core scattering is obtained from the diagram 
shown in Fig. 4. It evaluates to the S-wave t-matrix 


iTo(E) = exp {2i(Jo)ClD{E) , (20) 

with the Gamow-Sommerfeld factor C'^ = T{1 + ir])r[l — irf). 

In the absence of the Coulomb interaction, the t-matrix is usually ex¬ 
pressed in terms of effective range parameters. However, it is not possi¬ 
ble to separate the strong interaction from the Coulomb interaction in a 
model-independent way. Therefore one uses the so-called Coulomb-modihed 
effective range expansion (ERE) [20] to relate the phase shifts to redehned 
effective range parameters. The t-matrix is then written 


To(E) 


27r exp {2iao) 

mR A;C2(cot (5o — i) ’ 


( 21 ) 


where the total phase shift is given by (To-|-(5o- The S-wave Coulomb-modihed 
ERE is 


/cCAcot (5o - i) -f 2kcH{r]) = - — + -f ... , (22) 

Oq / 

where Oq and tq are the Coulomb-modihed scattering length and ehective 
range, respectively. One should note that the imaginary part of 2kcH{r]) 
exactly cancels —ikC^. 



Comparing Eqs. (20) and (21), order by order in the momentum k, we 
can express the scattering length and the effective range in terms of the 
low-energy coupling constants g and A 


1 

Uq 


27r 


5(2mR 


(A + S^'") , 


^^0 = - 




(23) 

(24) 


Equation (23) shows how the low-energy constant A absorbs the divergent 
part of the irreducible self-energy Furthermore, at LO we effectively 

take u to zero. The two parameters g and A are then not independent. 

The residue at the bound state, E = —B, of the full halo propagator 
defines the wavefunction renormalization, or the LSZ residue, and is therefore 
required for the calculation of bound-state observables. It is given by 




rd(D-i)i 

-1 

1 

dE 


“ z/ + S'( 

E=-B 


(25) 


In terms of the effective range, we write this as 


Gtt/cc 1 

9 ^™r iih, kc) - 3kcro ’ 


(26) 


using the matching condition in Eq. (24). In writing Eq. (26), we have defined 
the function 




6kl d 
ttir dE 


H{ri) 


E=-B 


(27) 


with the binding momentum 7 = The expression Eq. (26) is valid 

at next-to-leading order (NLO), that is it includes the effective-range correc¬ 
tion. At leading order (LO) the wavefunction renormalization is given by the 
simpler expression 




LO 


E'(-B) 

6nkc 


r^|i7(7,/cc) 


(28) 

(29) 


that is Eq. (26) with tq set to zero. Note that the factor of g^ in the de¬ 
nominator of the wavefunction renormalization will always cancel with a 
corresponding factor in an unrenormalized matrix element. 
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If the ratio between the binding momentum and Coulomb momentum 
'y/kc is small we can expand the function i/( 7 , kc) according to 

»(T.*:c) = l-|ir + ^ + ... . (30) 

Thus, for systems where the separation 7 /cc is fulhlled, we can use 
H{'y, kc) —)■ 1 in all expressions. (See Ref. [15] for a similar expansion of 

HM-) 

In this paper, we will £x the LSZ residue ^ to a calculated or measured 
asymptotic normalization coefficient (ANC), A, or to experimental radiative 
capture data. The ANC is defined as the coefficient in the asymptotic bound 
state wavefunction 

wi{r) = AhT_i^_)+i/ 2 ( 27 r) , (31) 

where W is the Whittaker-W function. The LSZ residue is related to the 
ANC according to 

^ ^ gW^lni + 

At LO, the wavefunction renormalization is determined solely by 7 and kc, 
as can be seen in Eq. (29), and as such this could in principle be used to 
predict the LO ANC as 


Alo = 


6kr 


//( 7 , kQ 


-r(i + kch) 


(33) 


At NLO we are left with one undetermined parameter, tq, in the LSZ 
residue (26). Thus, at orders beyond LO we can use an ANC as input to fix 
Z and to predict the effective-range parameters. We will therefore define our 
NLO wavefunction renormalization in terms of the matching to the ANC, 
Eq. (32). The ratio to the LO residue is then obtained as 


Z H{^,kc)A^ 

Zlo Qkc[T{l + kc/'y)f 


(34) 


Note that the matching (32) is valid at any order in the power counting and 
as such there is no EFT error due to the non-inclusion of higher-order contact 
interactions. 
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For a given one-proton separation energy, the ANC determines the Coulomb- 
modihed effective range 


^^0 


1 

3kc 


H{7,kc) 


6kcT{l + kc/'yf 

A 2 


(35) 


combining Eqs. (26) and (32). The Coulomb-modihed scattering length is 
then obtained from the pole position of the t-matrix (21), that is 


" 4kcH{-tkch)+'y^ro- ' ' 

These S-wave effective-range parameter predictions are accurate up to cor¬ 
rections of the shape parameter in the ERE. 


3. Charge form factor 

Information on the electromagnetic structure of an object can be obtained 
through elastic electron scattering. The elastic scattering amplitude can then 
be expressed in terms of the electric and magnetic form factors. Here we will 
focus on the electric charge form factor of an S-wave halo nucleus. The charge 
form factor, Fq, and charge radius, rc, are dehned by 


FciQ) - , n(P I-^emIp) 

e(Zc -h 1) 

(37) 

=i-^g2^... 

6 

(38) 


where p (p') is the incoming (outgoing) momentum of the halo field, Q = 
p' — p is the momentum transfer, and is the electromagnetic current. We 
evaluate this observable in the Breit frame, where no energy is transferred 
such that the photon four-momentum is (0, Q). 

In this paper, we calculate the charge form factor to NLO, by evaluating 
the relevant diagrams to this order. At LO there are two loop diagrams, 
rLo(Q), and at NLO a constant tree-level diagram, Tnlo; enters. We derive 
and evaluate these diagrams below. The charge form factor is then given by 
the sum of diagrams 

FciQ) = [rLo(Q) + rNLo(g) + ...]• (39) 
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Figure 5: The diagrams for the charge form factor at LO. The notation is as in Fig. 1. 
3.1. Leading order 

At leading order, we have to consider the diagrams shown in Fig. 5. 
There, the photon couples to the single-particle lines only, that is through 
an operator ip^A^ipk of dimension 5. We choose incoming and outgoing to¬ 
tal four-momenta as {E,—Q/2) and {E,Q/2), respectively. The resulting 
amplitude in momentum space is given by 

^rLo(Q) =t9^eZ,j 

X iSi{E - ho, -k3 + Q/2) txiks - /Q/2, k 2 - /Q/2, -B) 

X iSolho, k2) iSi{E — ho, —k2 -|- Q/2) 

X iSi{E - ho, -k 2 - Q/2) ix(k 2 -h /Q/2, ki -h /Q/2, -B) 

X iSoiho, ki) iSi{E - ho, -ki - Q/2) 

+ [(/ ~^ 1 ~ f), (Zc —)■ 1) , {So F-)- Ai)] , (40) 

where / = mo/Mtot is a mass ratio introduced for convenience. The last row 
makes sure that the photon also couples to the proton. Note that the total 
energy E is given by = —B + Q‘^/{SMtot), using energy conservation. 

We evaluate the energy integrals in Eq. (40) using the residue theorem, 
noting that the poles are located at ho = —i£ + k^/(2mo) for n = 1,2,3. 
The result can be written using two-body propagators, dehned in Eq. (3), as 

„-P f c / 

^rLo(Q) '^9 J (27r)^ *^tot( 

X x(k3, k2 - /Q/2, -B) S,h-B, k2 - /Q/2) 

X Stot{-B, k2 + /Q/2) x(k2 + /Q/2, ki, -B) 

X Atot(-5,ki) 

+ [(/ “^ 1 ~ f), {Zc —t 1)] , (41) 
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and this expression can be simplified using the Coulomb Green’s function in 
Eq. (6), to replace the two-body propagators ^tot and the four-point function 
y. This leads to 


irLo(Q) = -i9^eZ,J {k3|Go(-B)|k2-/Q/2) 

X (k 2 + /Q/2|Gc(-B)|ki) 

+ [(/ ~^ 1 “ /)) i^c -t 1)] • (42) 


By performing a Fourier transform on each of the momentum-space bras 
and kets, we arrive at the coordinate-space integral 

/ (T 3^ 

X (0|G'c(-5)|r) exp(i/Q-r) (r|G'c(-5)|0) 

+ [(/ 1 — /), {Zc —)■ 1)] . (43) 


This integral is much more convenient to use than the rather involved mo¬ 
mentum-space integral in Eq. (42). In the integral (43) the diagram in Fig. 5 
is also visualized better. It consists of two Coulomb Green’s functions, that 
propagate the fields from separation zero to r and back from separation r 
to zero, respectively, and the current operator in between the propagators. 
Since the Coulomb Green’s functions have one end at zero separation, only 
the S-wave part of these will contribute. Therefore we expand the Coulomb 
Green’s function in partial waves 

OO 

(r'|Go(£)|r) = J](2; + l)Gg(£;r',r)P,(i'''t) , (44) 

1=0 


with the bound-state partial-wave projected Coulomb Green’s function 


. ^ r 7 Fi{ri, p') [iFi{p, p) + Gi{p, p)] 

271 


(45) 


We can now write the Green’s function as (see Appendix A for details) 


(0|Go(-B)|r) =G<°>(-B;0,r) 

^ _ mRr(l + kch) hE_fcc/^q/2(27r) 

27r r 


(46) 
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As such, the integral (43) can be written as 


ffLoW) = - + kchf I 

+ [(/ ~^ 1 ~ f): {Zc —t 1)] • 


jo{fQr)W_kc/'r,i/2{‘^irY 

( 47 ) 


This integral can be evaluated numerically in a straightforward way. The LO 
charge form factor can now be calculated through 


FciQ) 


LO 


1 

e(-^c + 1 ) 


^LOrLo(Q) , 


(48) 


and the LO charge radius is given in terms of the loop-integral rLo(Q) and 
the wavefunction renormalization ^lO) according to 


r 


2 

c 


LO 


3Zlo d^ 

e(Zc + 1) dg2 


rLo(g) 


Q=0 


(49) 


We will now show that the charge form factor is normalized correctly to 
1 at g = 0. Starting from the coordinate-space integral (43) at g = 0 and 
the spectral representation of the Coulomb Green’s function 


(0|Gc(i?)|r) 


d^p ^p(0)^p(r) 

(27r)^ E — p2/(2mR) -|- ie 


(50) 


we hnd that 


rLo(0) = -/e(Z, + l) / dV|(0|Gc(-5)|r)f 


— — g^e{Zc -f 1) 


d^p ^p( 0 )^;( 0 ) 

(27r)3 {E — p2/(2mR) -|- ieY 


=e(Z, + l)S'(-5) . 


(51) 


In the first step above, a Dirac delta was used, coming from the integration 
over the orthonormal Coulomb wavefunctions 

j dV ^7p(r)^/’pi(r) = (27r)^(5(^)(p - p) . (52) 

The correct normalization of the LO charge form factor now follows, com¬ 
bining Eqs. (48), (51) and (28). 
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Figure 6: The diagram for the charge form factor at NLO. The notation is as in Fig. 1. 


3.2. Next-to-leading order 

The contributions that enter at NLO are that the full LSZ residue (25), 
or (32), is to be used and that the NLO tree-level diagram, given in Fig. 6, 
enters, through the operator d^A^d of dimension 6. This diagram is simply 
given by the Feynman rule for an Aq photon coupling to the dicluster held 


*rNLO = + 1) . 


(53) 


Note that the NLO correction diagram (53) is independent of the momentum 
transfer Q. These contributions come with an effective-range correction. 

The NLO charge form factor is given by the sum of the diagrams up to 
NLO, according to Eq. (39). It is clear that the normalization of the charge 
form factor at Q = 0 is still correct, using Eq. (51) together with the NLO 
LSZ residue (25), the NLO diagram (53) and the formula (39). The charge 
radius is then given by the order part of the LO loop-integral (47) together 
with the full wavefunction renormalization (32). The resulting NLO charge 
radius result is 


3Z 


d^ 


rr. = 


+ 1 ) dQ^ 


rL0W)ln., 


Q=0 


(54) 


It is evident that the NLO charge radius squared is a factor 


Z Hh,kc)A^ 

Zlo 6A;cF(1 + kchy ' 


(34) 


larger (or smaller) than the LO result. 

At higher orders there are three types of corrections. Firstly, there are 
local short-range operators 'ipl [V^Aq — (9o(V ■ A)] of dimension 7, which 
enter with hnite-size contributions of the core and proton helds. Secondly, 
there is a local short-range operator [V^Aq — i9o(V ■ A)] d, of dimension 8, 
that comes in with an undetermined short-range parameter. Thirdly, there 
are photon couplings due to the minimal substitution of derivatives in higher- 
order contact interactions, with the leading N^LO contribution coming from 
the shape-parameter in the ERE. 
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4. Radiative capture 


In this section we consider radiative capture of a proton into a halo state. 
By including range corrections, we are able to compute the cross section for 
this process to high orders in the Halo EFT expansion. We use quantum 
numbers relevant for the capture process ^®0(p, 7 )^^F*, which require that 
the incoming particle pair has relative angular momentum / > 1. Specihcally, 
we will consider the El capture through an incoming P-wave. 

The differential cross section for radiative capture of non-relativistic par¬ 
ticles is 


da 

dH 


Stt^p 




(55) 


where oo is the energy of the outgoing photon, p is the relative momentum 
of the incoming particle pair and e* are the photon polarization vectors. The 
vector amplitude A is for the capture process with a vector photon Ai being 
emitted. Note that we are working in Coulomb gauge, where the relation 


£j ■ Q = 0 , 


(56) 


for a real photon with momentum Q, is fulhlled. 

We present our results in terms of the astrophysical S-factor, which is 
dehned as 

S{E) = E exp {2717])atotiE) , (57) 

with the incoming center-of-mass energy E and the total cross section atot- 
From the Lagrangian (1), we can write down three classes of diagrams 
for the radiative capture process. Two of these are identically zero since the 
incoming particle pair is in a relative S-wave. For completeness we will show 
this explicitly below. The remaining diagram has contributions from partial 
waves / > 1, but since the P-wave dominates at low energies we neglect all 
other partial waves. 


4-1. Leading order 

At leading order, we only consider diagrams where the photon couples 
to one of the single-particle lines, that is through an operator 
of dimension 5. The two radiative capture diagrams of interest are shown 
in Fig. 7. We will write them using the pure Coulomb t-matrix Tq [21]. 
The four-point function y dehned in Fig. 2 is directly proportional to the 
Coulomb Green’s function Gc and receives contributions from the Coulomb 
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Figure 7: The diagrams for radiative capture at LO. The notation is as in Fig. 1. 


t-matrix Tc as well as from the free propagation of the core-proton system 
(cf. Fig. 2). In the center-of-mass of mass of the halo and the radiated photon, 
the momentum space diagrams from Fig. 7 are given by 

d^k2 


iA = - ig^/Zi^o j 


(27r) 


f5'o(/i;2o, k2) iSi{E - k2o - w, -k2 - Q) 


X 


*x(k2 + /Q, P + /Q, -B) iSi (pV(2 mi) - w, -p - Q) 


X {-ieZc{-p)/mi) 

/ d^fci 

*x(k2 + /Q, ki -h /Q, -B) iSi{E - kw - u, -ki - Q) 

X (-teZc(-ki)/mi) 

X iSi{E - kio, -ki) iSo{kio, ki) (-fTc(ki, p)) 

— [(/—)■ 1 — /) , {Zc —)■ 1) , (mo F-)- mi) , (S'o F-)- S'!)] , (58) 


where we have separated the contributions from the Coulomb t-matrix Tq 
and the free propagation of the core-proton system in the initial state. The 
total energy flowing through the diagram is E = —B 4- a; -|- o;^/(2Mtot) 
and the factor \/Z\^q is from the wavefunction renormalization of the halo. 
Performing the kno residue integrals and introducing a convenient Dirac delta, 
the integral (58) can be written as 


iA = — ig\/Z] 


LO' 


-eZcf f d^kid^k2 




(27r) 


Stot{—B, k2 -|- /Q) x(k2 -|- /Q, ki -|- /Q, —B) 


X 


x^tot(-5,ki + /Q) ki 

(5^(p - ki) + [E- k^/(2mR)]Tc(ki, p) 

[(/ 1 - /) , {Zc -)■ 1 )] . 


(59) 
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Then, using the Lippmann-Schwinger relation [21] 


-ippiki) = S^{ki-p) + [E-kl/{2mn)] ^Tc(ki,p) (60) 

and replacing the four-point function x with Eq. (6) the integral (59) can be 
expressed in a simpler fashion; 

iA = - ig\fz ^—^ ^ ^ (k2|G'c(-5)|ki + /Q) ki ^/>p(ki) 

“ [(/ “^ 1 ~ /) ) (^c 1 )] (61) 

Performing Fourier transforms and using ki exp (iki • r 2 ) = —iW 2 exp (iki ■ r 2 ) 
we can write the integral (61) as 

iA = - [ dVidV2 Gc\-B-,0,pri) 'ipp{r2) 

rriR J 

X exp {-ifQ ■ ri) V2(5®(r2 - ri) 

- [(/ ^ 1 - /) , (Z, ^ 1)] . (62) 


The resulting integral can now be expressed as 


iA =g\/Z] 


LO' 


- eZJ 

rriR 


d^r G'[?^(—5; 0,pr) exp (—i/car cos6*) (V'^p(r)) 


[(/ -t 1 - /) , {Zc -)■ 1 )] . 


(63) 


Summing over all polarizations and doing the angular integration, the am¬ 
plitude squared in Eq. (63) evaluates to 


E i‘. ■ -^ 1 ' 


■ Q( A 1 • ,x47r^eZc/exp(iai) 

V is sm c^(cos 0 -f sm <p) - 

m-Rp 


X I dr G'[?^(-5;0,p)jo(/c^r)^ [rFi(/ccMpr)] 


d 


“ [(/ “^ 1 “ /) ) i^C -t 1 )] 


(64) 


This integral can be solved numerically, using the S-wave Coulomb Green’s 
function given in Eq. (A.7) and the regular Coulomb wave function dehned 
in Eq. (9). 
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Figure 8: The first of the two vanishing classes of capture diagrams at LO. The notation 
is as in Fig. 1. 



Figure 9: The second of the two vanishing classes of capture diagrams at LO. The notation 
is as in Fig. 1. 


The remaining two classes of diagrams to consider at LO evaluate to zero 
for our choice of reference frame and kinematics. This can be understood 
from the fact that they correspond to an incoming S-wave. The diagram 
shown in Fig. 8 evaluates to zero in the zero-momentum frame, since the 
photon coupling is proportional to the ingoing momentum p' = 0 of the 
dicluster field: 




X 


. ^z/e(Zc + 1) /w . X 

i\JZ -—-p iD[E, p 

jLltot 

/ Mu 

iSo{ko, k) iSi{E - ko, -k) iTc(k, p) 


=0 


(65) 


The momentum-space amplitude from the next diagrams, which are shown 
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in Fig. 9, is given by 


j iSoih,, ka) iS,{E -u- ho, -ka - Q) 

X ix(k3 + /Q, k2 + /Q, -B) iSoiho, k2) iSi{E - u - k2o, -k2 - Q) 
X (-ieZc(-k2)/mi) iSi{E - ho,-'^ 2 ) ix0^2,'^3,-B) 

X iSoiho, ki) iSi{E - ho, -ki) {-ig) iD{E, 0) {-ig) 

/ d^fc 

i5'o(/i;o,k) iS'i(E -/cq,- k) iTc(k,p) 

+ [(/—)■ 1 — /) , {Zc —)■ 1) , (S'o <H- S'!)] . (66) 

Doing the same simplifications as before, it can be shown that this amplitude 
is parallel to the photon momentum: 

= — igy/Z [ d^r Gh—B] 0, p) exp (—i/Q ■ r) 

J 

X [VG'c(-5;p,0)] {-ig)iD{E,d){-ig)'4)hd) 

+ [(/ “^ 1 “ /) ) {^c —t 1)] 

= -tg^^D{E,d)%hf-^ [ dV G'c(-i?;0,p) 

J 

X ^(2/ + l)i^ji{fujr)Pi{-Q ■ f) r-^Gc{-B; p, 0) 
i L 

+ [(/ “^ 1 ~ /) ) (^c —t 1)] 

=Qg^^/ZDiE,0)^/JhO)^ [ dV Gc(-i?; 0 , p)ji(/a;r) 

mji J 

X ■^Gc{-B-,p,0) 

+ [(/ ~^ 1 “ /) ) i^c —t 1)] (67) 

In the first step the exponential function is expanded in spherical Bessel 
functions and in the second step all partial waves but I = 1 integrates to 
zero. Thus, using Eq. (56), 

Ci ■ = 0 (68) 

and the diagram does therefore not contribute. 
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Figure 10: The capture diagram entering at N'^LO. 


4-2. Next-to-leading order 

The higher-order ERE parameters appear with V -|-ieQA operators that, 
in principle, can give contributions to the radiative-capture amplitude. How¬ 
ever, the diagrams with these higher-order ERE operators are diagrams with 
initial-wave scattering due to the strong force. Since we only have included 
the S-wave part of the strong interaction these initial-wave scattering dia¬ 
grams are identically zero, which can be understood from the fact that the 
El capture process changes the angular momentum by one. If we were to 
include also the P-wave interaction explicitly in the held theory, then the 
effective range, the shape parameter, and so on, would contribute through 
diagrams with initial P-wave scattering. As the held theory is constructed 
in this paper, however, the physics of these diagrams is implicitly included 
in local short-range operators with growing powers of the photon energy uj. 
Such an operator is explicitly discussed below. 

Consequently, there are no additional capture diagrams to consider at 
NLO. The only contribution at NLO is due to the change in the wavefunction 
renormalization. This leads to a constant factor 

V^^LO 


larger (or smaller) result than the LO result. It is important to note that if 
ro ~ H{'j, kc)/{Skc), or equivalently if the ANC A is very large, then the 
NLO correction will be large, too. We discuss this in more detail in Sec. 5. 

The next correction enters at N^LO and as such this calculation is valid 
up to N^LO. 


H{'y, kc) - Skcro 
H{7,kc)A^ 
6kcT{l -h kc/'ff 


(69) 

(70) 
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J^.3. N^LO 

At N^LO a short-range El operator appears. The interaction term is 
given by 

£(n.lo) = [r>,vlc“ct + d\, - v,a„) (vA.^o,.) + hx., 

(71) 

and it has dimension 9. This short-range interaction is simply a contact 
vertex, where the incoming proton-core pair is in a relative P-wave. 

The operator (71) gives rise to the capture diagram in Fig. 10. The 
tree-level amplitude is given by 

B = D^^^^v^caexp (1 + ri^)C^ . (72) 

The derivation of the amplitude (72) involves the evaluation of the P-wave 
integral J j|^k'^p(k) = — exp (i(Ti)py^ (1 -|- where Ui is the pure 

Coulomb phase shift in the P-wave (see App. B of Ref. [14] for details on 
P-wave integrals). The symbol has been introduced as a compact no¬ 

tation for the total constant of proportionality. 

The next relevant operator 

i:(N.LO) = [h/7’C.“Ct + 49o (d,>A, - V,X1„) (^i^.^„,,)+hx., 

(73) 

of dimension 11, enters at N®LO and as such this calculation is valid up to 
N^LO. 

5. Fine tuning and S-wave proton halos 

Along the neutron drip-line there exist several neutron halo states. These 
states are characterized by an unnaturally large neutron-core scattering length, 
which brings the state very close to threshold. However, proton halos are 
much more rare. In the S-wave case this can be understood by considering 
the Coulomb repulsion between the valence proton and the core. For a proton 
halo state to exist, we need the attractive strong force to be almost cancelled 
by the Coulomb repulsion, resulting in a threshold state [22, 23]. This can¬ 
cellation can be seen within our formalism as an additional hne-tuning in 
the effective range. Note that for an S-wave proton halo state this means 
that the proton-core scattering length needs to be unnaturally large and that 
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Figure 11: The dependence of the LO charge radius on kd^. The solid black line is the 
LO Halo EFT result, the blue circle denotes the ^^F* system, and the dashed red line is the 
asymptotic 1/fcc behavior. The inset shows the low-energy region in a semi-logarithmic 
scale, illustrating the hypothetical neutron halo limit (74). The curve was generated using 
a binding momentum 7 = 13.6 MeV. 


the effective range must be fine-tuned to cancel the Coulomb repulsion. The 
existence of proton halos is therefore doubly suppressed by the need for two 
fine-tunings. 

Proton halo systems contain the second scale fee that will depend on the 
charge of the degrees of freedom in the system. Within our framework, kc 
is a parameter that is independent of the Coulomb-modified effective range 
parameters. When /cc /7 3> 1, we will speak of the extreme Coulomb regime. 
In this regime the two particles tend to be close together, since otherwise the 
system is ripped apart by the Coulomb repulsion. This limit is in practice 
realized for the system, where /cc/t = 3.8. In Fig. 11 the LO charge 
radius is shown as a function of the Sommerfeld parameter /cc/t, where we 
have used the binding momentum for ^^F*, 7 = 13.6 MeV. The blue circle 
is the parameter point corresponding to the physical ^^F* system. It is clear 
that this system is almost in the extreme Coulomb regime. Note that the 
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resulting LO charge radius is very small for a strong Coulomb repulsion, 
since it has an asymptotic 1/kc behavior. At the far left, where kc 7 , the 
system mimics that of a neutron halo, with the only difference being that 
the photon also can couple through minimal substitution to the nucleon held. 
The limiting value for A;c /7 —)■ 0 is therefore given by [7] 


lim Tq = 


1 ZJ^ + (1 - ff 

Zc + 1 272 


(74) 


In the standard power counting for systems with large S-wave scattering 
length the effective range enters at NLO. The hierarchy of scales in this case 
is 7 ,/cc < 1/ro l/i?core, where Rcore IS the length scale set by the core. 
However, the discussion above implies that we can have /cc^o ~ 1 instead 
of kc In the zero-range limit, the inverse Coulomb momentum sets 

the scale for the LO charge radius and the effective range contributions will 
therefore be numerically large since the LSZ-factor for S-wave proton halo 
nuclei with the effective range included behaves as 


Z oc 


1 

1 - 3kcro 


(75) 


It appears to be fine tuned to the pole position with ro ~ l/(3A;c)- In the 
case of the system the effective-range correction results in a factor 3.6- 
3.8 larger charge radius. This hierarchy-of-scales problem can be solved by 
hxing the bound state pole position of the t-matrix at leading order and the 
ANC at NLO. This procedure ensures that corrections beyond NLO scale 
naturally again. 


6. S-factor for ^®0(p, 7 )^^F* 

The excited state of is due to an S-wave interaction between the 
valence proton and the ^^0(0^) core. As such, we describe this proton halo 
state using the Halo EFT formalism presented in this paper. 

As was discussed in Sec. 5, the Coulomb momentum is larger than the 
binding momentum for ^^F*, that is kc 3 > 7 - Consequently, the effective- 
range correction is needed for reliable predictions. The ANC for this sys¬ 
tem has been extracted by Huang et al. [25], using a single-particle model 
fit of radiative capture data, as Ai = 77.21 fm“^'^^ and experimentally by 
Gagliardi et al. [26], using the transfer reaction ^®0(^He,(i)^'^F, as A 2 = 
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E (keV) 


Figure 12: Energy dependent S-factor for 7 )^^F* fitted to experimental data by 

Morlock et al. [24]. The error bands correspond to the model error from omitted terms at 
higher orders. See text for details. 


(80.6 ± 4.2) Comparing these ANCs with the LO result (33), Apo = 

21.4 makes it clear that the effective range must be very close to the 

pole position l/{3kc), such that the LSZ residue becomes large. 

In Fig. 12 we present our Halo EFT results for the radiative proton cap¬ 
ture reaction ^®0(p, 7 )^’^F* together with data by Morlock et al. [24]. The 
green dashed line is the Halo EFT result valid up to N^LO. The single free 
parameter of this model is the effective range, and this was fitted to the ex¬ 
perimental data by Morlock et al. [24] by minimizing an objective function 
defined as 



of-(a) - C'r'‘ T 

Si ) 


(76) 


where Of^ and denote the theoretical and experimental values of ob¬ 
servable Oi in the pool of fit data M, and the total uncertainty 6i determines 
the weight of the residual. The theoretical values, and therefore the resid¬ 
uals, depend on the vector of fit parameters a. The total uncertainty is 
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the squared sum of experimental and theoretical errors, Sf = 5exp,i + 

The truncation of the Halo EFT expansion allows to estimate the model 
error. At N^LO the omitted terms should scale as {p/khi)^, where p is the 
incoming momentum and /chi the breakdown scale of the theory. At this 
order, we therefore assign a theoretical error (5th,i = C'eft {Pi/khi)^- Taking 
these energy-dependent model errors into account will allow us to include 
relatively high energy data in the £t. We use data up to a center-of-mass 
energy of 2.3 MeV and determine the amplitnde, Ceft, of the model error 
by the statistical gniding principle that the total per degree of freedom 
should be unity. Note that the Morlock data has a normalization error of 
10%, which we snbtracted in quadratnre from the total experimental error 
dnring the fitting procedure. This normalization error was added back after 
the fit had been performed. The constant Ceft is expected to be of natnral 
size and is determined iteratively snch that the degree of freedom is 

minimized to nnity. This systematical theory error at N^LO is shown as a 
green band in Fig. 12, with Ceft = 6.9 and the breakdown scale given by 
khi = 76 MeV. Althongh somewhat large, this valne of Ceft is still consistent 
with onr power connting estimate. 

Similarly, the red line is given by a fit of the N^LO model to the same 
data set where the systematic theory error estimate scales as (p/Zchi)®- At 
this order we find Ceft = 1-9 with k^i = 76 MeV. It is clearly seen that 
the result converges with increasing order of the EFT and that the S-factor 
valne at threshold is stable. From these fits, we extract a threshold S-factor 

^ _ f (9.9 ± 0.1 (stat) ± 1.0 (norm)) keV b , N^LO , . 

I (l0.4 ± 0.1 (stat) ± 1.0 (norm)) keV b , N^LO 

with the 1% error due to the EFT fit (mainly statistical error) and the 10% 
error from the nncertainty in the absolnte normalization of the experimental 
data. These results give the ANC 

A — I (77-4 ± 0.2 (stat) ± 3.8 (norm)) fm“^A ^ N^LO ^ 

( (79.3 ± 0.2 (stat) ± 3.9 (norm)) fm“^A ^ N^LO ’ 

which is consistent with the ANCs of Huang et al. and Gagliardi et ai. 

The charge radius of the ^^F* can now be obtained by using an extracted 
ANC. Using the ^®0-proton ANC extracted from the N^LO radiative proton 
capture fit, the resulting NLO charge radius is given by 

rc,NLO = (2.20 ± 0.04 (EFT) ± 0.11 (ANC)) fm . (79) 
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The NLO EFT error in Eq. (79) was estimated from the EFT expansion 
parameter squared, (q/Zchi)^, using a breakdown scale /chi = 76 MeV, which is 
of the same order as the inverse core radius l/i?core ~ 73 MeV. The dominant 
error in Eq. (79) is from the normalization error of the Morlock data, through 
the extracted ANC. However, there could also be additional EFT errors for 
this result due to the non-inclusion of the operators that are responsible for 
the hnite-size contributions of the constituents [27]. Alternatively, we can 
interpret our result as the radius relative to the core. 

7. Conclusions 

We have calculated the charge radius and radiative capture cross section 
for proton halo nuclei interacting through a large S-wave scattering length. 
Specihcally, we have included higher-order effective-range corrections and 
shown consequent good agreement with experimental data. Our description 
of proton-capture on agrees very well with the data and leads to a new 
way of extracting the low-energy S-factor and the corresponding ANC with 
error estimates that take the intrinsic effective theory error into account. For 
the excited state of ^^F we found a large correction to all observables at NLO 
and have explained this as a result of the complicated interplay between the 
range of the strong interaction and the scale set by the Coulomb interaction. 
We showed that constraining LO and NLO to proton bound state pole po¬ 
sition and ANC will warrant that higher order corrections scale naturally. 
Our results highlight that Halo EFT provides a powerful tool to analyze ex¬ 
perimental data and predict observables with theoretical error estimates. A 
further advantage of the analysis of such systems with an EFT framework 
is the proper identihcation of two-body contributions to the electromagnetic 
current operator. The physics of such terms should in principle also appear 
in a cluster model description although there are some contributions that 
are not generated by gauging the strong interaction terms and are gauge 
invariant by themselves. The EFT power counting also provides a reliable 
estimate for the relative importance of these terms. 

Our analysis can also be used to obtain the low-energy phaseshifts for 
nucleon-nucleus elastic scattering that can then be compared to correspond¬ 
ing ab initio calculations that have become possible over the last years [28]. 

It is a very interesting question, whether a two-proton halo could in prin¬ 
ciple support an effective three-body bound state and how the spectrum 
of this system compares to two-neutron halo nuclei bound by resonant S- 
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wave interactions. This question is also relevant for studies of two-proton 
radioactivity [29, 30]. Furthermore, it seems worthwhile to investigate the 
description of low-lying resonances in proton-|-core systems. 
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Appendix A. Partial wave projected Coulomb Green’s function 

In this appendix, we discuss the construction and partial wave projection 
in the S-wave channel. 

It is convenient to use the Coulomb Green’s function in a non-integral 
form and below we present such a form for the bound state Green’s function. 
This can be done by doing a partial-wave projection and forming the Green’s 
function as a product between two independent Coulomb wavefunctions, sat¬ 
isfying one boundary condition each, in accordance with the definition of 
Green’s function. For the r = 0 boundary condition we must use the regular 
Coulomb wave function Fl and to satisfy the condition for a bound state at 


r = oo we need to form the combination 



iFi + Gl . 

(A.l) 

This 

can be seen from the asymptotics 



Fl(? 7 , p) sin (p - L7r/2 - p log (2p) -h ai) 

(A.2) 

and 

Gl{v, p) cos (p - LtiI2 - plog (2p) -h uz,). 

(A.3) 



using that for a bound state p = i 7 r, with 7 > 0 , where the only combination 
that yields only an exp (—yr) dependence is the combination given in (A.l). 
Therefore, the partial wave projected Coulomb Green’s function is 




m^pFL{p,p') [iFL{p,p) +Gl{p,p)] 
271 p'p 


(A.4) 


where the normalization is given by the Wronskian of the Coulomb wave- 
functions and the Coulomb-Schrodinger equation. 

For the S-wave bound state observables we consider, only propagation 
with a Coulomb Green’s function down to zero separation is needed, that is 
i?;0,r). Using that 


iFi( 7 ,p)+G'L( 7 ,p) = y I exp + ip)jfU_i^,L+i/ 2 (- 2 ip) , 

(A.5) 

which was shown in [19], and the limit 

lim = exp (— 7rp/2)A/r(l -|- ip)r(l — ip) (A. 6) 

p^O p 


we may write Eq. (A.4) for S-waves as 
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P 


(A.7) 
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